Creation of solitons and vortices by Bragg reflection of Bose-Einstein condensates in 

an optical lattice 
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We study the dynamics of Bose-Einstein condensates in an optical lattice and harmonic trap. 
The condensates are set in motion by displacing the trap and initially follow simple semiclassical 
paths, shaped by the lowest energy band. Above a critical displacement, the condensate undergoes 
Bragg reflection. For high atom densities, the first Bragg reflection generates a train of solitons and 
vortices, which destabilize the condensate and trigger explosive expansion. At lower densities, soliton 
and vortex formation requires multiple Bragg reflections, and damps the center- of- mass motion. 
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Optical lattices (OLs) provide unprecedented control 
of transport through the energy bands of periodic quan- 
tum systems. This has led to beautiful experimental 
demonstrations of Bloch oscillations [l]] and quantized 
Wannier- Stark ladders Q for non-interacting ultra-cold 
alkali atoms. There is also great interest in understand- 
ing the behavior of Bose-Einstein condensates (BECs), 
formed from interacting alkali atoms, in OLs [p|-|l5|]. Pre- 
dictions that accelerated condensates will perform 
Bloch oscillations, whose turning points at the top of 
the energy band correspond to successive Bragg reflec- 
tions, have been confirmed in experiments f7|,|l2[ on 
BECs with equilibrium peak densities no % 10 14 cm 
For no > 10 14 cm -3 , more complex motion has been 
observed , which cannot be explained by Bragg re- 

flection or analogous semiclassical models of energy band 
transport. Previous numerical studies of condensate dy- 
namics in OLs have used the one-dimensional (ID) Gross- 
Pitaevskii equation [||J|Jl^]. They provide invaluable in- 
sights for understanding the mean center-of-mass motion 
of the condensate, but have not related this motion to 
changes in the internal structure of the condensate, in 
particular dynamical excitations such as solitons and vor- 
tices. Producing such excitations in a controlled way 
requires state-of-the-art experimental techniques, which 
involve manipulating the condensate phase and/or den- 
sity profile |l5| - |l9[ , rotating the confining trap [20,21], 
moving a laser beam through the atom cloud |l3|,22,23], 
or tuning the inter-atomic interactions p4| . 

In this Letter, we show that Bragg reflection provides 
a new mechanism for generating solitons and vortices 
in BECs. Moreover, these excitations can have a dra- 
matic effect on the evolution of the atom cloud. At the 
first Bragg reflection, the condensate wavefunction is a 
standing wave with nodes at each maximum in the OL 
potential. At each node, the condensate phase changes 
abruptly by tt. Bragg reflection therefore imprints atom 
density and phase profiles similar to those used to gener- 
ate solitons in experiment [p~5|— p~9[ . The effect of this im- 



printing on the condensate dynamics depends critically 
on the atom density. For condensates with no > 10 14 
cm -3 , realized in recent experiments pHJl^| , it leads to 
the self-assembly of a chain of stationary solitons, which 
decay rapidly into vortex anti- vortex pairs. Strong inter- 
actions between the vortices destabilize the atom cloud, 
causing it to explode and fragment. For no < 10 14 
cm -3 , the standing wave formed at the first Bragg re- 
flection produces no dynamical excitations. But subse- 
quent Bragg reflections do generate solitons and vortices, 
which damp the center-of-mass motion. The dissipation 
and instability processes that we identify could play a 
key role in the complex dynamics recently observed for 
high-density condensates in OLs p| , |l2| . 

We consider condensates formed from Na 87 Rb atoms 
in a ID OL and a three-dimensional harmonic trap. 
Figure 1(a) shows the potential energy profile of the 
OL, Vol(x) = ^o sin 2 (7nr/<f), whose depth Vb = 23 
peV and period d = 397.5 nm are taken from experi- 
ment [|ll| . The trap frequency for confinement along the 
z— direction is high enough for the BEC dynamics to re- 
duce to two-dimensional (2D) motion with potential en- 
ergy V(x, y) = VoL(x) + \m{u'lx 2j cUyy 2 ) 1 where m is the 
mass of a single atom and uo x1 uo y are frequencies of the 
harmonic trap. For most of the results presented here, 
u x = 2tt x 50 rad s -1 , u y = 2tt x 35 rad s _1 , Na = 10 4 
, and no = 0.43 x 10 14 cm -3 (System A). We also con- 
sider a second set of parameters, uo x = 2tt x 8.7 rads -1 , 
Uy = 2tt x 90 rads -1 , and iV A =3x 10 5 (System B), 
corresponding to recent experiments |ll[. In this case, 
no is sufficiently high (~ 1.5 x 10 14 cm -3 ) for the stand- 
ing wave formed during Bragg reflection to have a par- 
ticularly dramatic effect on the BEC. For both sets of 
parameters, uo x is small enough to ensure that the har- 
monic potential energy variation across each OL period 
is much less than the width Y = 10 peV of the lowest en- 
ergy band (Fig. 1(a)). Consequently, the harmonic trap 
only weakly perturbs the band structure [|5| . 

We determine the density profile of the condensate by 
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FIG. 1. (a) Solid curve: potential energy profile of the OL. 
Grey rectangle: energy range of the lowest energy band. Dot- 
ted curve: x— dependence of the harmonic potential energy, 
(b) Initial density profile of System A, along y = 0. Inset: 
density profile and x— dependence of harmonic potential en- 
ergy (dotted) immediately after trap displacement. 

using the Crank-Nicolson method || to solve the time- 
dependent Gross-Pitaevskii equation 
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where a — 5.4 nm is the s- wave scattering length [|26| , and 
^(x,y,t) is the condensate wavefunction at time £, nor- 
malized so that |\I/| 2 is the number of atoms per unit vol- 
ume |2?[ . The equilibrium density profile for the ground 
state of System A is shown in Fig. 1(b). 

At t — 0, we disrupt the equilibrium of the condensate 
by suddenly displacing the harmonic potential through 
a distance Ax along the x— axis jnj]. Displacing the 
trap increases the initial potential energy of the BEC by 
AV ~ ^mujl(Ax) 2 (Fig. 1(b) inset). As the atoms start 
to move, this potential energy is converted into kinetic 
energy, which determines how far the BEC has acceler- 
ated up the lowest energy band |25[| . In order for the 
condensate to reach the top of the band and therefore 
undergo Bragg reflection, AV must be > T, which, for 
System A, requires Ax > (2T/mcj x r )^ = 15/im. 

We now consider the condensate dynamics obtained 
from Eq. (1) for Ax = 10/im, below the threshold for 
Bragg reflection. Figure 2(a) shows that the mean 
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FIG. 2. (a) (x) versus t for System A with Ax = 10 /im. 
Insets: grey-scale plots of density (black high) in x — y plane 
(axes inset) at t =10.7 ms (left) and 21.3 ms (right), (b) 
Grey-scale plot: \f(k x ,t)\ 2 (white = 0, black high) for System 
A with Ax = 10/xm. Open circles: points on corresponding 
semiclassical trajectory k x {t). (c) Solid curve: (x) versus t for 
System A with Ax = 25/xm. Arrow marks first turning point. 
Dashed curve: corresponding semiclassical orbit x(t). (d) As 
(b), but for Ax = 25/xm. (e) (x) versus t for System B with 
Ax = 150/xm. Insets: grey-scale plots of density (black high) 
in x — y plane at t = 13.3 ms (top) and 18 ms (middle). Lower 
inset shows enlargement of boxed region in upper inset, (f) 
As (b), but for System B with Ax = 150/xm and omitting the 
k x (t) curve which, for t < 13 ms, is indistinguishable from the 
narrow Fourier distribution. 

(center-of-mass) position of the condensate, (x), under- 
goes simple periodic motion, bounded by the harmonic 
trap. The internal structure of the BEC is unaffected by 
this motion, being the same at t = 10.7 ms (Fig. 2(a) 
left inset) and t = 21.3 ms (Fig. 2(a) right inset) as 
at t = 0. To determine how the condensate moves in 
reciprocal space, we calculate the Fourier transform of 
Sfr(x,0,t). The Fourier power, \f(k x , t)| 2 , corresponding 
to wavenumber k Xl remains narrow and changes period- 
ically as t increases (grey-scale plot in Fig. 2(b)). Since 
the condensate's internal structure does not change with 
t when Ax = 10/im, the form of \f(k x , t)\ 2 and the corre- 
sponding (x) versus t curve (Fig. 2(a)) can be understood 
by considering the motion of a single point particle in 
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FIG. 3. Grey-scale plots of density (white = 0, black high) 
in x — y plane (axes inset) for System A with Ax = 25 /xm and 
t = ms (a), 7.5 ms (b), 10.7 ms (c). Plots are symmetrical 
under reflection about y = 0. Vertical dotted lines indicate 
x = in each case. Horizontal bar shows scale. Cross in (b) 
marks center of a soliton. Region within dashed box in (c) is 
shown in Fig. 5(a). 

the lowest energy band. The single particle trajectories 
x(t) and k x (t) in real and reciprocal space are deter- 
mined by the semiclassical equations of motion dx/dt = 
h~ 1 dE(k x )/dk x and dk x /dt = h~ 1 F x p8||, where E(k x ) is 
the energy- wave number dispersion relation for the band, 
and F x = —mcj^x is the harmonic restoring force along 
the x— direction. In Fig. 2(b), the Fourier power (grey- 
scale plot) is concentrated along the single-particle k x (t) 
curve (white circles). The corresponding real-space tra- 
jectory, x(t), is indistinguishable from the plot of (x) ver- 
sus t shown in Fig. 2(a). 

When Ax is increased to 25/im, above the threshold for 
Bragg reflection, the mean x position of the condensate, 
determined from Eq. (1), performs damped periodic mo- 
tion (solid curve in Fig. 2(c)). In Fig. 2(d), the grey-scale 
plot of \f(k x ,i)\ shows that the condensate' s mean k x 
value increases approximately linearly with t and reaches 
the Brillouin zone boundary at 2.6 ms. At this time, 
the condensate undergoes Bragg reflection, which pro- 
duces the first (arrowed) turning point in Fig. 2(c). 
The quantum calculations of (x) and \f(k x ,t)\ 2 deviate 
rapidly from the corresponding semiclassical trajectories, 
x(t) and k x (t), shown respectively by the dashed curve 
and open circles in Figs. 2(c) and (d). For t > 7.5 ms, 
the oscillations in (x) are damped and multiple peaks ap- 
pear in \f(k x , t)| 2 , which spreads through the Brillouin 
zone. This deviation from single-particle behavior indi- 
cates that the BEC's center-of-mass motion is strongly 
affected by changes in its internal structure. Key stages 
in the evolution of the 2D density profile are shown in 
Fig. 3. As t increases from (Fig. 3(a)), the density 
minima deepen and fall to zero at the first Bragg reflec- 
tion, which we now analyze in detail. 

The lower curves in Figs. 4(a), (b), and (c) show 
|^(x,0,t)| 2 just before (t = 2 ms), at (t = 2.6 ms), and 
just after (t = 3 ms) the first Bragg reflection. The up- 
per curves show the wavefunction phase, modulo 
2tt. Just before reflection (Fig. 4(a)), the density near 
the center of the BEC has a minimum value of ~ 10 13 
cm -3 , which is approximately half that at t = (Fig. 
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FIG. 4. Lower curves: density profiles along y = for the 
condensate in System A with Ax = 25 /xm and t = 2 ms (a), 
2.6 ms (b), 3 ms (c). Upper curves show <j>{x) modulo 27T, 
with vertical scale indicated by bars of length tt. 

1(b)). The local velocity along the x— direction, v x = 
(h/m)d(f)/dx, is > throughout the condensate. At the 
point of Bragg reflection (Fig. 4(b)), the density min- 
ima fall to zero at each peak in Vol(x). At each zero, 
(j) changes abruptly by tt (upper curve in Fig. 4(b)). 
Away from the discontinuities, d(j)/dx ~ 0, indicating 
that the BEC is at rest. This variation of density and 
phase demonstrates that a standing wave forms at the 
point of Bragg reflection. In recent experiments, laser 
illumination was used to produce similar individual den- 
sity minima and/or tt phase shifts [|l5]-|l9|, which sub- 
sequently evolved into dark solitons. By analogy, the 

standing wave might be expected to generate a chain of 

_ i 

stationary solitons, each of width w « (iranM) 2 , where 
tim is the local mean atom density J15[ . At the first 
Bragg reflection, um — 3 x 10 13 cm~^near the center 
of the condensate in System A (see Fig. 4(b)), and so 
w ~ 3.5<i [jnj. Since w is much larger than the width 
(~ d) of the density minima in the standing wave, the 
first Bragg reflection does not produce solitons in Sys- 
tem A. Instead, after reflection, the density minima rise 
away from zero and d(j)/dx becomes negative for all x, as 
the condensate starts to move from right to left. 

We now consider the condensate motion for t > 5 
ms. Figure 3(b) shows the density profile at the sec- 
ond Bragg reflection, when t = 7.5 ms. Again a standing 
wave is formed, which creates nodal lines in the den- 
sity profile (white stripes in Fig. 3(b)) at each maxi- 
mum in Vol(x). However, in contrast to the first Bragg 
reflection, the standing wave now disrupts the internal 
structure of the condensate sufficiently to allow a soli- 
ton, marked by the cross in Fig. 3(b), to form across 
several OL periods. Subsequent Bragg reflections gener- 
ate more solitons, which have a pronounced effect on the 
condensate's internal structure and center-of-mass mo- 
tion. To illustrate this, Fig. 3(c) shows the density profile 
at t = 10.7 ms. For x > 0, there are three extended soli- 
tons (white crescent shapes), whose wavefronts have been 
curved by refraction originating from the non-uniform 
density [O] . This refraction is the precursor of snake in- 
stability |13|], which causes the solitons to decay into two 
vortices with opposite circulation, like those enclosed by 
arrows in Fig. 3(c). At the center of each vortex, |^| 2 
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FIG. 5. (a) Grey-scale plot of density for System A within 
dashed box in Fig. 3(c) (white = 0, black high). Arrows show 
direction of circulation around vortices, (b) Grey-scale plot 
of <fi (white — 0, black — 2tt). Ax — 25 /im and t — 10.7 ms. 

= 0. The vortices can be seen more clearly in Figs. 
5(a) and (b), which show respectively enlargements of 
the density profile within the dashed box in Fig. 3(c), 
and the corresponding phase. The soliton represented by 
the white crescent in Fig. 5(a) appears as a dark island 
towards the right of Fig. 5(b). At the left-hand edge of 
this island, the phase changes abruptly from 7r/4 (light 
grey) to 37r/4 (dark grey). Around the two vortices, the 
phase (Fig. 5(b)) changes continuously from (white) to 
2tt (black), indicating quantized circulation in the direc- 
tion of the arrows. Vortex formation is the main cause of 
damping in the center-of-mass motion (Fig. 2(c)). A cru- 
cial aspect of this damping mechanism is that the soliton 
formation and subsequent vortex shedding occur when 
the condensate is almost at rest. It is therefore funda- 
mentally different to the phonon emission process used 
to interpret experiments on System B Q], which occurs 
when v x exceeds a critical value of ~ 5 mm s -1 . It is 
also unrelated to the damping found in ID simulations of 
Bloch-oscillating condensates |^| , which cannot include 
the effects of vortex formation. 

We now relate our calculations to the experiments on 
System B [[ll|. Figures 2(e) and (f) show the time evo- 
lution of (x) and \ f(k x ,t)\ 2 for this system, after a large 
trap displacement of 150/im. As in System A, the first 
Bragg reflection generates a density node and associated 
7r phase shift at each maximum in Vol(x). But since um 
is much larger in System B ( ~ 3.3 x 10 14 cm -3 near 
the BEC center at the first Bragg reflection), w ~ d. 
Since w is so closely matched to the width of the density 
minima in the standing wave, Bragg reflection causes the 
self-assembly of ~ 10 stationary solitons, which form a 
chain across the central third of the condensate. Figure 
2 (e) shows the compact cigar-shaped density profile of 
the condensate just after the first Bragg reflection (up- 
per inset). The region within the box is shown enlarged 
in the lower inset of Fig. 2(e), which reveals three sta- 
tionary solitons (extended white areas). The time taken 
for solitons to develop following phase imprinting is pro- 
portional to the distance / over which the phase changes 
by 7r p5|. In a standing wave, / ~ 0, and so the first 



Bragg reflection and the formation of the soliton chain 
occur almost instantaneously. The solitons decay rapidly 
into two chains of vortex ant i- vortex pairs, which form 
a complex interacting system. The interactions create 
a massive internal strain, which causes the BEC to ex- 
plode laterally (perpendicular to the x-axis) , resulting in 
the diffuse and fragmented atom density profile shown in 
the middle inset of Fig. 2(e) p0[ . The explosion has a 
dramatic effect on the /^—distribution of the atoms (Fig. 
2(f)), which is initially extremely narrow but, at the point 
of Bragg reflection, spreads through the whole Brillouin 
zone. This could account for the broad momentum distri- 
bution observed when high-density condensates undergo 
Bragg reflection Q . 

In summary, we have investigated how Bragg reflection 
affects the internal structure and center-of-mass-motion 
of condensates accelerating through an OL. When the 
atom density is high enough to ensure that w < d, 
the density zeros and tt phase shifts imprinted by the 
first Bragg reflection generate a train of stationary soli- 
tons, which decay rapidly into vortex ant i- vortex pairs. 
Strong interactions between the ensemble of vortices have 
a catastrophic effect on the condensate, causing it to un- 
dergo explosive expansion. This dynamical regime is a 
unique feature of condensates in an OL and should be 
experimentally accessible in existing systems fll| , |l2| • For 
lower atom densities, soliton formation requires multiple 
Bragg reflections. The subsequent decay of the solitons 
into vortex anti- vortex pairs provides a new dissipation 
mechanism, which could contribute to the damping of the 
center-of-mass oscillations observed in experiment [[111. 
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